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drag  and  ground  friction.  Predicted  values  of  maximum  velocity, 
displacement  at  maximum  velocity,  and  total  displacement  were 
tabulated  for  1224  exposure  conditions.-  These  conditions  included 
both  air  and  surface  bursts  with  yields 'from  1 T to  100  MT  and  peak 

overpre.saures  from  2 to  97  psi. 

^)®iological  criteria  were  presented  which  indicate  that  per- 
sonnel subjected  to  decelerative  tumbling  over  open  terrain  can 
tolerate  much  higher  velocities  than  personnel  impacting  a non- 
yielding, flat  surface  at  normal  incidence..  These  criteria  were 
used  in  preparing  range-yield  plots  for  1-,  2.5-.  5-,  and  50- 
percent  probabilities  of  serious  injury  to  prone  personnel  in  the 
vicinity  of  a surface  hurst. 

thods  for  extending  the  presented  results  to  other  exposure 
conditions  were  discussed 
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Thi?5  work  unit,  Dccelorative  Tumbling  Blast  Effects 
Personnel.  Contract  No.  DNA  00j_-74-C-0120,  was  part  of  the 
subt?,..-Dk,  Biological  Effects  of  Blast  Environments. 

This  study  was  under  the  direction  of  R.  K.  Jones. 
The  technical  assistance  of  W.  Hicks,  K.  Saunders,  A.  Shaw, 
and  W.  S.  Jackson  is  acknowledged.  T.  Minagawa  performed 
the  technical  photography  and  prepared  the  illustrative  ma~ 
terial.  The  report  was  typed  and  compiled  by  B.  Martinez. 

This  research  was  conducted  according  to  the  prin- 
ciples enunciated  in  the  "Guide  for  Laboratory  Animal  Fa- 
cilities and  Care,"  prepared  by  the  National  Academy  of 
Sciences,  National  Research  Council. 
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INTRODUCTION 


One  of  the  major  hazards  to  personnel  In  the  vicinity 
of  an  explosion  is  associated  with  whole-body  displacement 
induced  by  blast  winds.  The  purpose  of  this  study  was  to 
predict  the  probability  of  impact  injury  from  displacement 
as  a function  of  yield  and  ground  range.  Predictions  were 
made  for  personnel  in  different  orientations  in  open  terrain 
and  near  structural  complexes.  A mathematical  model  (Ref- 
erence 1)  was  used  to  calculate  the  time-displacement  history 
of  personne]  from  considerations  of  aerodynamic  drag  and 
ground  friction.  The  model  had  been  verified  previously 
using  data  from  anthropomorphic  dummies  translated  over  open 
terrain  by  high-explosive  (HE)  and  nuclear  detonations 
(References  2 and  3). 

Two  types  of  terrain  were  of  interest  because  of  dif- 
ferences in  the  translational  velocities  required  to  produce 
injuries.  In  open  terrain,  personnel  undergo  decelerative 
rumbling,  in  which  case  the  energy  of  motion  is  dissipated 
slowly  through  a series  of  glancing  impacts,  and  relatively 
high  velocities  can  be  tolerated.  In  contrast,  personnel 
translating  near  a structural  complex  may  impact  rigid  struc- 
tures at  approximately  normal  inc idence , and  the  energy  of 
motion  is  dissipated  much  more  rapidly  with  a corresponding 
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increase  in  the  likelihood  of  injuries.  The  biological 
criteria  used  in  the  present  study  were  derive!  from  experi- 
ments (Reference  4)  by  airblast  from  a shockti.be.  Injuries 
commonly  incurred  during  either  type  of  deceleration  included 
lacerations,  contusions,  fractures,  dislocations,  and  damage 
to  internal  organs. 


MATHEMATICAL  MODEL 

The  mathematical  model  used  to  calculate  blast- 
induced  translation  of  personnel  (Reference  1)  v.as  adapted 
to  an  electronic  digital  computer.  The  computer  program  is 
quite  general  in  that  it  takes  into  account  variable  ground 
friction,  aerodynamic-drag  characteristics,  and  blast-wave 
parameters . 

Ground  Friction 

The  average  deceleration,  due  primarily  to  frictional 
effects  with  the  ground,  of  a 165-lb  man  undergoing  decelera- 
tive  tumbling  has  been  estimated  in  Reference  1 to  be  ap- 
pro.ximateJy  8.9029  ft/soc^,  where  Vj  i,s  the  in- 

stantaneous translational  velocity  in  ft/sec.  This 
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relationship  was  derived  from  translational  data  obtained 
with  guinea  pigs,  rabbits,  dogs,  and  a goat  by  scaling  to 
an  animal  the  size  of  man.  Previous  calculations  have 
indicated  that  the  ground- f ri ct ion  formula  also  applies 
reasonably  well  to  full-scale  anthropomorphic  dummies  (Ref- 
erence 3 ) . 

Acceleration  Coefficient 

In  order  to  calculate  the  aerodynamic  drag  forces  on 
a man  it  is  necessary  to  know,  at  each  instant,  the  accel- 
eration coefficient,  a , defined  as  the  projected  area  pre- 
sented to  the  wind  multiplied  by  the  drag  coefficient  and 
divided  by  the  body  mass  (Reference  5).  The  a varies  as  the 
orientation  of  a man  changes  during  translation.  a value  of 
approximately  0.03  ft^/lb  has  been  estimated  to  be  the  aver- 
age a for  a random  orientation  (i.e.,  considering  all  possi- 
ble orientations  to  be  equally  likely)  of  an  unbent  man. 

The  a for  a prone  man  end-on  to  the  wind  is  0.0063 
ft^/lb  (Reference  6).  Data  obtained  on  500-ton  HE  tests 
(Reference  3)  suggest  that,  for  near-ideal  conditions,  asyme- 
tries  will  cause  a dummy  which  was  initially  prone  and  end- 
on  to  the  flow  to  turn  side-on  to  the  wind  after  about  13  ft 
of  displacement . The  details  of  the  subsequent  displacement 
were  not  obtained,  but  the  assumption  of  a random  orientation 
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seems  reasonable.  An  emperic.al  fit  to  the  data  yielded  the 
following  relationships; 

a = 0.0063  + (Di/31)4  for  0 ^Df<12.2 
a = 0.03  for  12.2  4 Df 

where  a is  expressed  in  ft^/lb  and  Df  is  the  instantaneous 
displacement  in  feet . 

The  a for  a prone  man  in  a random  orientation  i.s 
0.014  ft2/lb  (Reference  6).  This  value  is  smaller  than  the 
average  a for  a man  considering  all  orientations  because  a 
prone  man  cannot  be  face-  or  back-on  to  the  wind.  The  data 
mentioned  in  Reference  3 suggest  that,  on  the  average,  a 
dummy  which  was  initially  prone  in  a random  orientation  will 
turn  side-on  to  the  wind  after  about  6 ft  of  displacement. 

The  following  relationships  were  used  to  approximate  the 
acceleration  coefficient  during  the  displacement  oi  a man 
initially  prone  in  a random  orientation: 

a = 0.014  + (Df/47)2  for  0 4Di<5.9 
a.  = 0.03  for  5.9  4 Oi 

The  a for  a prone  man  side-on  to  the  blast  is  0.022 
ft^/lb.  The  data  indicate  that  dummies  initially  in  this 
orientation  will  start  to  roll  within  a very  short  displacement 
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(^1  ft).  The  average  a for  a man  rolling  with  his  long  axis 
perpendicular  to  the  wind  is  0.038.  It  Is  not  known  after 
what  displacement  the  orientation  would  typically  become 
more  random  in  nature.  Because  of  these  uncertainties,  an 
ol  of  0.03  ft‘^/lb  was  assumed  throughout  the  displacement  of 
a man  initially  prone  and  side-on  to  the  wind. 

The  only  calculations  made  for  a standing  man  were 
for  the  worst  case,  i.e.,  for  a man  initially  face-  or  back-on 
to  the  wind.  Data  obtained  with  dummies  indicate  that  these 
two  initial  orientations  will  result  In  approximately  the 
same  displacement  (Reference  3).  It  should  be  noted  that 
orientation  probably  influences  the  displacement  of  a stand- 
ing man  less  than  in  the  case  of  a prone  man,  in  that  the  a cf  a 
standing  man  can  range  from  only  0.022  to  0.052,  whereas  the 
a of  a prone  man  can  range  from  0.0063  to  0.022  depending  on 
the  oriontation. 

At  overpressures  below  2.7  psi , the  standing  dummies 
toppled  over  with  their  feet  having  undergone  almost  no  dis- 
placement. At  overpressures  above  5.3  psi,  the  legs  were 
blown  out  from  under  the  dumndes  such  that  they  became  air- 
borne almost  immediately.  Typically  a standing  dummy  rotated 
such  that  it  first  impacted  the  ground  on  its  head  after  a 
displacement  of  from  10  to  15  ft.  Rather  than  attempting  to 
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adapt  the  mathematical  model  to  closely  approximate  the  con- 
plexities  introduced  by  the  sometimes  airborne  displacements, 
a simpler  approach  was  considered.  Calculations  were  made 
assuming  that  the  ct  remained  constant  at  0.052  and  that  the 
ground  friction  was  operative  throughout  the  displacement 
(Reference  3).  The  predictions  agree  well  with  experimental 
displacements  obtained  over  a limited  range  of  yield  (50  T 
to  38  KT)  and  overpressure  (5  to  50  psi).  The  good  agree- 
ment is  probably  due  in  part  to  two  compensating  errors: 

(1)  assuming  the  a remains  constant  at  0.052  should  result  in 
the  predicted  velocities  being  too  high  because  the  dummies 
rotated,  thus  reducing  the  a;  (2)  assuming  the  ground  friction 
to  always  be  operative  should  result  in  the  predicted  veloci- 
ties being  too  low  because  the  dummies  were  airborne  during 
the  first  10  ft  of  displacement. 

Blast-Wave  Parameters 

The  blast-wave  parameters  used  to  predict  the  dis- 
placement of  personnel  were  derived  primarily  from  data  for 
nuclear  detonations  in  References  7 and  8.  Calculations  were 
made  at  selected  yields  and  ground  ranges  for  surface  bursts 
and  for  upper  and  lower  optimum  heights-of-burst  (HOB's). 

The  optimum  HOB  is  defined  as  the  one  which  tends  to  maximize, 
for  a given  ground  range,  both  the  peak  dynamic  pressure  and 
the  impulse  of  the  positive  dynamic  pressure.  Thus  the  optimum 
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IIOB  would  also  tend  to  maximize  the  translational  velocity 
of  personnel  for  a given  ground  range.  The  airblast  data 
indicate  that  at  intermediate  ground  ranges  there  are  two 
local  optimum  HOB's,  with  the  upper  optimum  HOB  giving  the 
larger  dynamic-pressure  levels  and  impulses  for  scaled 
ground  ranges  greater  than  approximately  770  ft/(KT)^^^, 
and  the  lower  optimum  HOB  giving  the  larger  dynamic-pressure 
levels  and  impulses  for  smaller  scaled  ground  ranges. 

Table  A-1  list.s  the  pertinent  blast-wave  parameters. 
Note  that  values  are  given  for  both  upper  and  lower  optimum 
HOB'S  for  scaled  ground  ranges  between  400  and  950  ft/(KT)l/3. 
For  surface  bursts  the  peak  overpressures  given  in  the  table 
range  from  2 to  97  psi  corresponding  to  peak  dynamic  pres- 
sures between  0.093  and  140  psi. 

Double-exponential  equations,  of  the  types  used  in 
Reference  8,  were  fitted  to  the  overpressure  and  dynamic- 
pressure  waves  to  facilitate  their  being  entered  into  the 
computer  program.  The  Rankine-Hugon iot  relations  were  assumed 
to  apply  at  the  shock  front,  and  changes  within  the  shock  wave 
were  assumed  to  be  adiabatic.  It  was  thus  possible  to  use  the 
overpressure  to  calculate  the  air  density,  which  in  turn  could 
be  combined  with  the  dynamic  pressure  to  calculate  the  wind 
speed,  all  of  these  quantities  varying  with  time. 
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In  calculating  a displacement,  the  blast-wave  param- 
eters were  evaluated  at  fixed  time  intervals,  the  size  of 
which  depended  on  the  yield  of  interest.  For  yields  of  1 T, 

1 KT , 1 MT , and  100  MT  the  time  intervals  were  5 psec,  50 
ysec , 0.5  msec,  and  2 msec,  respectively.  These  values  also 
represent  the  time-step  sizes  used  in  solving  the  differential 
equations  of  motion  by  Runge-Kutta  methods. 

CALCULATED  VELOCITIES  AND  DISPLACEMENTS 

Each  of  the  34  runs  listed  in  Table  A-1  was  analyzed 
in  connection  with  nine  yields  (ranging  from  1 T to  100  MT 
by  factors  of  10)  and  four  initial  orientations  (three  prone 
and  one  standing).  Thus  a total  of  1224  displacements  were 
predicted  using  the  computer.  Figure  A-1  shows  a typical 
velocity-displacement  history  for  a prone  n.rn  in  a random 
orientation  exposed  to  an  overpressure  level  of  27.4  psi 
from  a 10-KT  surface  burst.  As  was  the  usual  case,  the 
maximum  velocity,  V,  was  achieved  within  a displacement,  S, 
which  was  short  compared  to  the  total  displacement,  D,  Like- 
wise, the  time  to  peak  velocity  was  short  compared  to  the 
total  time  of  displacement. 

The  predicted  V,  S,  and  D values  for  all  of  the  com- 
puter runs  are  listed  in  Tables  A-2  through  A-5.  For  a given 
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yield  and  initial  orientation,  the  V and  D values  always  in- 
creased with  the  peak  overpressure.  The  S values  also  in- 
creased with  the  peak  overpressure  except  for  some  slight 
reversals  at  the  higher  overpressures.  These  reversals  may 
be  an  anomaly  caused  by  the  finite  size  of  the  time  step 
used  in  the  calculations.  Even  if  this  is  the  case,  the  ef- 
fect on  the  ""omputed  V and  D values  would  be  minimal.  The 
predicted  maximum  velocities  in  Tables  A-2  through  A-5  for 
personnel  exposed  to  a surface  burst  have  been  plotted  in 
Figures  A-2  through  A-5. 

INJURY  PROBABILITY  VS  RANGE  AND  YIELD 

Table  A-6  indicates  the  probability  of  serious  injury 
(fracture  or  ruptured  internal  organ)  as  a function  of  ve- 
locity for  personnel  undergoing  either  decelerative  tumbling 
or  impact  at  normal  incidence  against  a nonyielding,  flat 
surface.  The  velocitie.s  are  in  agreement  with  estimates  of 
man's  tolerance  to  impact  in  Reference  9. 

Figures  A-2  through  A-4  were  used  to  determine  the 
scaled  ground  ranges  corresponding  to  the  velocities  in  Table 
A-6.  Ground  range  vs  nuclear  yield  was  plotted  in  Figures 
1 and  2 for  1- , 2.5-,  and  5-percent  probabilities  of  injury 
to  prone  personnel  and  in  Figures  3 and  4 for  50-percent 


13 


► 5%  PROBABILITY  OF  IMPACT  INJURIES 


GROUND  RANGE,  R.  ft 

- to  5-Percent  Probability  of  Impact  Injuries  to  Prone  Personnel 
n Ocen  Terrain.  The  right-hand  edge,  center,  and  left-hand  edge 
f each  curve  correspond  to  probabilities  of  1,  2.5,  and  5 percen 
espectively.  Computed  for  a surface  burst. 


50%  PROBABILITY  OF  IMPACT  INJURIES 
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10  MT 


Figure  4.  50-Percent  Probability  of  Impact  Injuries  tc  Prone  Personnel  Near 
Structures.  Computed  for  a surface  burst. 


probabilitiy  of  injury.  Implicit  in  Figures  2 and  4 is  the 
assumption  that  personnel  near  structures  impacted,  at  nor- 
mal incidence,  a nonyielding,  flat  surface  at  the  instant  of 
ma.ximum  velocity. 

Some  of  the  curves  in  Figures  1 through  4 do  not  ex- 
tend to  1 T because  no  attempt  wa.s  made  to  extrapolate  to 
yields  less  than  those  that  could  be  read  directly  from  the 
curves  in  Figures  A-2  through  A-4 . However,  the  missing 
portions  of  the  curves  correspond  to  peak  overpressures  of 
greater  than  55  psi,  such  that  the  primary-blast  hazard 
would  probably  override  the  displacement  hazard. 

From  Figures  1 through  4 it  can  be  seen  that,  ior  a 
given  injury  probability  and  yield,  the  ground  range  is  very 
dependent  on  both  initial  orientation  and  typo  of  terrain. 
For  example,  for  a 50-percent  probability  of  injury  from  a 
10-KT  surface  burst,  the  ground  ranges  for  end-on  and  side- 
on  prone  personnel  are  approximately  850  and  1500  ft  in  open 
terrain  and  1400  and  2500  ft  near  structures.  In  contrast, 
varying  the  injury  probability  from  1 to  5 percent  has  a 
relatively  small  influence  on  ground  range.  For  example, 
for  a 10-KT  surface  burst  in  open  terrain,  the  ground  ranges 
for  1-  and  5-percent  probabilities  of  injury  are  approxi- 
mately 1200  and  1100  ft  for  end -on  prone  personnel  and  2000 
and  1900  ft  for  side-on  prone  personnel . 
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Four  methods  of  extending  the  possible  applications 
of  the  predictions  presented  in  this  paper  are  worth  noting: 
First,  velocities  read  from  Table  A-5  or  Fig- 
ure A-5  could  be  used  to  estimate  injury  probability 
vs  range  and  yield  for  standing  personnel. 

Second,  it  would  be  possible  to  predict  tlie 
approximate  displacement  of  personnel  subjected  to 
a near-classical  dynamic-pressure  wave  generated  by 
an  unspecified  source  by  assuming  that  (for  a given 
orientation)  V,  S,  and  D are  each  a function  of  only 
the  peak  dynamic  pressure  and  the  impulse  of  the 
positive  dynamic  pressure.  Further,  as  the  duration 
of  the  wave  decreases,  the  significance  of  the  wave 
shape  and  peak  dynamic  pressure  decreases,  such  that 
V and  D become  functions  of  only  the  impulse  of  the 
positive  dynamic  pressure,  and  S approaches  zero. 

All  of  these  functions  could  be  derived  from  the  data 
in  Tables  .A-1  through  A-5. 

Third,  for  a HOB  lower  than  those  listed  in  Table 
A-1,  it  is  reasonable  to  use  logarithmic  interpolation 
to  estimate  translational  velocities  and  displacements. 
For  example,  for  a HOB  of  200  ft/(KT)^/^,  one  could 
interpolate  between  values  for  a surface  burst  and  a 
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lower  optimum  HOB  for  ground  ranges  between  400  and 
950  ft/(KT)^/^,  and  interpolate  between  values  for 
a surface  burst  and  an  upper  optimum  HOB  for  ground 
ranges  between  1100  and  2GOO  ft/(KT)^'^^. 

Finally,  the  predictions  may  be  applied  directly 
to  an  HE  detonation  by  considering  a nuclear  yield 
twice  as  large  as  the  HE  yield.  Thus,  for  500  tons 
of  high  explosives,  one  would  read  1 KT  on  the  nu- 
clear yield  scale.s  in  the  figures.  This  procedure 
i-s  based  on  tlio  assumption  that  approximately  one- 
half  of  the  energy  released  in  a nuclear  detonation 
goes  into  airblast  . 
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Figure  A-4 . Maximum  Translational  Velocity  of  Frone  Personnel  Side-On  to 
the  Wind.  Computed  for  a surface  burst. 


Figure  A-5.  Maximum  Translational  V'elocity  of  Standing  Personnel  Front 
Back-On  to  the  Wind.  Computi  d for  a surface  burst. 


TABLE  A-1 


Uax  Dynamic  Proflsur<',  psi 
Gr^mnd  Rangp,  ft 

Duration  of  Positive  Overpressure,  wee 
Duration  <jf  Positive  Dynamic  Preswure.  sec 


TABLE 


TABLE  A-4 
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